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ABSTRACT 

We examine volume-limited samples of galaxies drawn from the Sloan Digital Sky Survey to look 
for relations among internal and collective physical parameters of galaxies as faint as Mr = —17.5. 
The internal physical properties of interest include morphology, luminosity, color, color gradient, con- 
centration, size, velocity dispersion, equivalent width of Ha line, and axis ratio. Collective properties 
that we measure include the luminosity and velocity dispersion functions. We morphologically classify 
galaxies using the three dimensional parameter space of u — r color, g — i color gradient, and concentra- 
tion index. All relations are inspected separately for early and late type galaxies. At fixed morphology 
and luminosity, we find that bright {Mr < —20) early type galaxies show very small dispersions in 
color, color gradient, concentration, size, and velocity dispersion. These dispersions increase at fainter 
magnitudes, where the fraction of blue star-forming early types increases. Late types show wider 
dispersions in all physical parameters compared to early types at the same luminosity. Concentration 
indices of early types are well correlated with velocity dispersion, but are insensitive to luminosity and 
color for bright galaxies, in particular. The slope of the Faber- Jackson relation (L c>c a"^) continuously 
changes from 7 = 4.6 ± 0.4 to 2.7 ± 0.2 when luminosity changes from Mr = —22 to —20. The size of 
early types is well-correlated with stellar velocity dispersion, tr, when a > 100 km s^^. We find that 
passive spiral galaxies are well separated from star-forming late type galaxies at Ha equivalent width 
of about 4. An interesting finding is that many physical parameters of galaxies manifest different 
behaviors across the absolute magnitude of about ± 1. The morphology fraction as a function of 
luminosity depends less sensitively on large-scale structure than the luminosity function (LP) does, 
and thus seems to be more universal. The effects of internal extinction in late type galaxies on the 
completeness of volume-limited samples and on the LP and morphology fraction are found to be very 
important. An important improvement of our analyses over most previous works is that the extinction 
effects are effectively reduced by excluding the inclined late type galaxies with axis ratios of b/a < 0.6. 
Subject headings: galaxies: general - galaxiesduminosity function, mass function - galaxies:formation 
- galaxies:fundamental parameters - galaxies:statistics 



1. INTRODUCTION 

With the large and homogeneous redshift surveys like 
the Sloan Digital Sky Survey (SDSS; York et al.2000) 
and Two Degree Pield Galaxy Redshift Survey (2dPGRS; 
Colless et al. 2001), our view of the local universe be- 
came extended out to a few hundred mega parsecs. Prom 
such dense redshift surveys, the small-scale distribution 
of galaxies is also being revealed in more detail. It is now 
possible to make very accurate measurements of various 
relations among the physical properties of galaxies and 
the local environment. This is a difficult task, which 
should be completed before theoretical modeling and in- 
terpretations are attempted. The study is complicated 
because there are many physical parameters involved and 
their mutual relations are non-trivially interrelated. A 
further complication is the fact that definition of envi- 
ronment depends on the type of objects chosen to trace 
the large scale-structure. 

One could divide the physical properties of galaxies 
into internal and collective ones. Examples of internal 
physical properties are color, luminosity, morphology, 
star formation rate, velocity dispersion, surface bright- 
ness, spectral type, size, and so on. Collective properties 
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can be the strength and statistical nature of spatial clus- 
tering, peculiar velocity field, luminosity function, veloc- 
ity dispersion distribution function, halo mass function, 
axis ratio distribution function and so on (see Pig. 1). On 
the other hand, the local environments can be spatial or 
temporal. Dependence on the latter means the evolution 
in cosmic time. "Spatial environment" can have various 
meanings. Traditionally discrete types, such as cluster, 
group, field, or void, are used to distinguish among dif- 
ferent environments of galaxies. Recently, the smoothed 
number density of neighboring galaxies was also used to 
represent the environment. One could also use other pa- 
rameters to define the environment of galaxies. Density 
gradient, shear fields, and large-scale peculiar velocity 
fields can be important parameters affecting the proper- 
ties of galaxies. The relations among all these physical 
observables and environment can give us information on 
galaxy formation and the background universe. Many 
authors in the past have investigated the relations be- 
tween environment and galaxy properties (Lewis et al. 
2002; Gomez et al. 2003; Goto et al. 2003a; Balogh et 
al. 2004a, 2004b; Hogg et al. 2004; Tanaka et al. 2004; 
Kuehn & Ryden 2005; Desroches et al. 2006). 

There has been much recent progress in studying the 
correlation between the distribution of galaxies and their 
basic properties, and the correlation among the proper- 
ties. Por example, the luminosity function of galaxies 
subdivided by various criteria has been measured. Blan- 
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Fig. 1. — Schematic diagram showing some key science for understanding the structure formation and the background universe. 



ton et al. (2003a) presented the luminosity functions of 

the SDSS galaxies defined by galaxy light profile shape. 
Madgwick et al. (2002) measured the luminosity function 
of 2dFGRS for different type of galaxies defined by their 
spectral properties. Nakamura et al. (2003) studied the 
dependence of the luminosity function on galaxy mor- 
phology and the correlation between galaxy morphology 
and other photometric properties such as color and con- 
centration index using a sample of SDSS galaxies with 
visually identified morphological types. Weinmann et al. 
(2006) performed an extensive study of the dependence 
of color, star formation, and morphology of galaxies on 
halo mass using galaxy groups selected in the SDSS DR2. 
Shen et al. (2003) examined the correlation between the 
size distribution of galaxies and other properties such 
as their luminosity, stellar mass and morphological type. 
Baldry et al. (2004) analyzed the bivariate distribution of 
the SDSS galaxies in the color versus absolute magnitude 
space. The fundamental plane and the color-magnitude- 
velocity dispersion relation for early type galaxies with 
SDSS have been investigated by Bernardi et al. (2003a, 
2005) Very recently Desroches et al. (2006) showed that 
the fundamental plane projections of elliptical galaxies 
depend on Imuinosity using DR4 of the SDSS. In partic- 
ular, they found that the radius-luminosity and Faber- 
Jackson relations are steeper at high luminosity, which 
we also found in our work using a larger SDSS sample 
with morphologically cleaner early type samples. Sheth 
et al. (2003) and Mitchell et al. (2005) measured the 
velocity dispersion distribution function of early type 



galaxies carefully selected from the SDSS samples us- 
ing selection criteria similar to those of Bernardi et al. 
(2003a). Alam & Ryden (2002) measured the axis ratio 
distributions for red and blue galaxies fitted by the de 
Vaucouleurs and exponential profiles. 

This paper focuses on understanding of the relation- 
ship among many physical properties of galaxies us- 
ing a set of volume-limited samples of the SDSS galax- 
ies, which are further divided into subsamples of con- 
stant morphology and absolute magnitude. This well- 
controlled experiment allows us to measure the relations 
among internal and collective properties as a function of 
morphology and luminosity. In our companion paper we 
have studied the dependence of various galaxy proper- 
ties on the local density environment (Park et al. 2006, 
hereafter Paper II). 

This paper is organized as follows. In section 2, we 
describe our catalog and morphology classification and 
the physical parameters that we use in Section 3 and 4 
to investigate the relations among the galaxy properties. 
We summarize our results in Section 5. 

2. OBSERVATIONAL DATA SET 

2.1. Sloan Digital Sky Survey 

The SDSS (York et al. 2000; Stoughton et al. 2002; 
Adelman-McCarthy et al. 2006) is a survey to explore 
the large-scale distribution of galaxies and quasars by 
using a dedicated 2.5m telescope at Apache Point Ob- 
servatory (Gunn et al. 2006). The photometric survey 
has imaged roughly tt sr of the northern Galactic cap in 
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five photometric bandpasses denoted by u, g, r, i, and 
z centered at 3551,4686,6165,7481, and 8931^1, respec- 
tively, by an imaging camera with 54 CCDs (Fukugita 
et al. 1996; Gunn et al. 1998). The limiting magni- 
tudes of photometry at a signal-to-noise ratio of 5 : 1 
are 22.0, 22.2, 22.2, 21.3, and 20.5 in the five bandpasses, 
respectively. The median width of the point-spread func- 
tion (PSF) is typically 1.4", and the photometric uncer- 
tainties are 2% rms (Abazajian et al. 2004). 

After image processing (Lupton et al. 2001; Stoughton 
et al. 2002; Pier et al. 2003) and calibration (Hogg et 
al. 2001; Smith ct al. 2002; Ivczic ct al. 2004; Tucker et 
al. 2006), targets arc selected for spectroscopic follow- 
up observation. The spectroscopic survey is planned to 
continue through 2008 as the Legacy Survey and to yield 
about 10^ galaxy spectra. The spectra are obtained by 
two dual fiber-fed CCD spectrographs. The spectral res- 
olution is A/AA ^ 1800, and the rms uncertainty in red- 
shift is ~ 30 km s~^. Because of the mechanical con- 
straint of using fibers, no two fibers can be placed closer 
than 55" on the same tile. Mainly due to this fiber col- 
lision constraint, incompleteness of spectroscopy survey 
reaches about 6% (Blanton et al. 2003b) in such a way 
that regions with high surface densities of galaxies be- 
come less prominent even after adaptive overlapping of 
multiple tiles. This angular variation of sampling density 
is accounted for in our analysis. 

The SDSS spectroscopy yields three major samples: 
the main galaxy sample, the luminous red galaxy sample 
(Eisenstein et al. 2001), and the quasar sample (Richards 
et al. 2002). The main galaxy sample is a magnitude- 
limited sample with an apparent Petrosian r-magnitude 
cut of TOr,iim ~ 17.77, which is the limiting magnitude for 
spectroscopy (Strauss ct al. 2002). It has a further cut 
in Petrosian half-light surface brightness fimo, fim = 24.5 
mag/ arcsec^. More details about the survey can be found 
on the SDSS web site. ^ 

In our study of galaxy properties, we use a subsam- 
ple of SDSS galaxies known as the New York University 
Value-Added Galaxy Catalog (NYU-VAGC; Blanton et 
al 2005). This sample is a subset of the recent SDSS 
Data Release 5. One of the products of the NYU-VAGC 
used here is a large-scale structure sample DR4plus (LSS- 
DR4plus). For local density estimation in the three- 
dimensional redshift space in our accompanying paper 
we use galaxies within the boundaries shown in Figure 1 
of Paper II, which improves the volume-to-surface area 
ratio of the survey. There are also three stripes in the 
southern Galactic cap observed by SDSS. Density esti- 
mation is difficult within these narrow stripes, so for con- 
sistency with the samples examined in the accompanying 
paper, we do not use them. The remaining survey region 
covers 4464 deg^. The primary sample of galaxies used 
here is a subset of the LSS-DR4plus sample referred to 
as "voidO", which is further selected to have apparent 
magnitudes in the range 14.5 < r < 17.6 and redshifts 
in the range 0.001 < z < 0.5. These cuts yield a sample 
of 312,338 galaxies. The roughly 6% of targeted galaxies 
that do not have a measured redshift due to fiber colli- 
sions are assigned the redshift of their nearest neighbor. 

Completeness of the SDSS is poor for bright galaxies 
with r < 14.5 because of both the spectroscopic selection 

^ http://www.sdss.org/dr5/ 



criteria (which excludes objects with large fiux within 
the 3" fiber aperture; the cut at r = 14.5 is an em- 
pirical approximation of the completeness limit caused 
by that cut) and the difiiculty of obtaining correct pho- 
tometry for objects with large angular size. For these 
reasons, analysis of SDSS galaxy samples have typically 
been limited to r > 14.5. This is unfortunate because 
it limits the range of luminosity that can be probed; us- 
ing the magnitude limits of the voidO sample the range 
of absolute magnitude is only about 3.1 at a given red- 
shift. To extend the range of magnitude, we attempt to 
complete the redshift sample by supplementing the cat- 
alog with redshifts of bright galaxies, first from SDSS 
itself, and then by matching to earlier redshift catalogs 
(the missing galaxies are quite bright and easily observed 
by earlier surveys). We employ another subset of the 
LSS-DR4plus sample referred as bvoid, which is made 
before redshift determination and includes galaxies with 
r < 14.5. Galaxies in bvoid with SDSS redshifts are 
added to the catalog. For those objects lacking red- 
shifts, we searched earlier catalogs. Out of 5227 bright 
objects with r < 14.5 in the bvoid sample and within 
our survey region, we assign redshifts to 3945 galaxies 
cross matched with the Updated Zwicky Catalog (UZC; 
Falco et al. 1999), 807 galaxies from the bvoid sample 
after redshift determination, 69 galaxies from the Cen- 
ter for Astrophysics (CfA) redshift catalog'' (ZCAT 2000 
Version), and 4 galaxies from the IRAS Point Source 
Catalog Redshift Survey (Saunders et al. 2000). We 
found 336 bright objects to be errors, such as outlying 
parts of a large galaxies that were deblended by the auto- 
mated photometric pipeline, blank fields, stars, etc. Out 
of the remaining 66 objects, 30 galaxies are included in 
our sample, using redshifts from the NASA Extragalac- 
tic Database (NED^). Only 36 galaxies remain in the 
bright bvoid sample with no measured redshift. In ad- 
dition to these bright bvoid galaxies, which have been 
selected based on their r-band magnitudes, we have in- 
cluded the redshifts of 340 UZC galaxies cross-matched 
with SDSS objects in NYU-VAGC, which were not in- 
cluded in the previous step. Out of 6295 UZC galaxies 
within our survey region, 242 are still not matched with 
SDSS objects in the NYU-VAGC. Ofthe missing UZC 
galaxies, 60% have redshifts smaller than our inner red- 
shift cut of z = 0.025, leaving only 99 galaxies missing 
at larger redshift. 

In total, we add 5195 bright (r < 14.5) galaxies to 
the voidO samples, which yields a final sample of 317533 
galaxies. Volume-limited samples derived from the re- 
sulting catalog have nearly constant comoving number 
density of galaxies along the radial direction at redshift 
z > 0.025. Thus, we treat our final samples as having no 
bright limit at redshift greater than z = 0.025. 

2.2. Definitions for the Volume-limited Samples 

In this study we use only volume-limited samples of 
galaxies selected by absolute magnitude and redshift lim- 
its. Figure 2 shows the definitions of our six subsam- 
ples in redshift-absolute magnitude space. The shallow- 
est subsample is labeled CM, which stands for 'color- 
magnitude'. We employ the CM sample to extend some 

* http;//www.cfa.harvard.edu/~huchra/zcat/zcom.htm 
^ http://nedwww.ipac.caltech.edu 
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Fig. 2. — Sample definitions of our six volume-limited SDSS samples in redshift vs. absolute magnitude space. 



of our analyses to fainter absolute magnitude. The 
deeper and thicker (along the line of sight) five samples, 
from Dl to D5, are the main samples for study of galaxy 
properties. The definitions for all samples are summa- 
rized in Table 1. 

The comoving distance and redshift limits of each 
volume-limited sample defined by an absolute magnitude 
limit are obtained by using the formula 

TOrjim - M^jim = 5log{r{l + z)) + 25 + Kiz) + E{z), (1) 

where K{z) is the mean X-correction, E{z) is the mean 
luminosity evolution correction, and r is the comoving 
distance corresponding to redshift z. We adopt a flat 
ACDM cosmology with density parameters JIa = 0.73 
and flm = 0.27 to convert redshift to comoving distance. 
To determine sample boundaries we use a polynomial fit 
to the mean i^-correction, 

K{z) = 3.0084(z-0.1)2-Kl.0543(z-0.1)-2.51og(l + 0.1). 

. .(2) 

We apply the mean luminosity evolution correction given 
by Tegmark et al. (2004), E{z) = 1.6(z - 0.1). The 
rest-frame absolute magnitudes of individual galaxies are 
computed in fixed bandpasses, shifted to z = 0.1, using 
Galactic reddening corrections (Schlegel et al. 1998) and 
X-corrections as described by Blanton et al. (2003c). 
This means that a galaxy at z = 0.1 has a X-correction 
of — 2.51og(l+0.1), independent of its spectral energy dis- 
tribution (SED). We have applied the same mean lumi- 
nosity evolution correction formula to individual galax- 
ies. 



The physical parameters we consider in this study are 
'^■^{u — r) color, absolute Petrosian magnitude in the r- 
band M^, morphology (see below), Petrosian radius, axis 
ratio, concentration index, color gradient in "'^(g — i) 
color, velocity dispersion, and equivalent width of the 
Ha line. 

To compute colors, we use extinction and /^-corrected 
model magnitudes. The superscript 0.1 means the rest- 
frame magnitude i^-corrected to the redshift of 0.1. All 
our magnitudes and colors follow this convention, and 
the superscript will subsequently be dropped. To mea- 
sure some of the physical parameters of galaxies, we 
retrieve the g- and i-band atlas images and basic pho- 
tometric parameters of all galaxies in our sample from 
the archive of photometric reductions conducted by the 
Princeton/NYU group. ^ To take into account flattening 
or inclination of galaxies, we use elliptical annuli in all 
parameter calculations, and the isophotal position angle 
and axis ratio in the i-band are used to define the ellip- 
tical annuli. The g — i color gradient is defined to be 
the color difference between the region with R < O.bRpd 
and the annulus with O.SiZpct < R < Rpot, where i?pct is 
the Petrosian radius. The (inverse) concentration index 
is defined by R^o/Rqq where i?5o and Rgo are the semi- 
major axis lengths of ellipses containing 50% and 90% of 
the Petrosian fiux in the i-band image, respectively. The 
g — i color gradient, concentration index, and isophotal 
axis ratios are corrected for the effects of seeing. To do 
so, we first generate a large set of Sersic model images 
with various scaling lengths, slopes, and axis ratios con- 
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TABLE 1 
Volume-limited Samples 



Name 


Absolute Magnitude 


Redshift 


Distance ^ 


GaIaxies(AfB'') 




CM 


-17.5 > Mr 


0.025 <z < 0.03484 


74.6 <R< 103.7 


11756 (3467) 


3.00 


Dl 


-18.0 > Mr 


0.025 <z < 0.04374 


74.6 <R< 129.9 


20288 (6256) 


3.41 


D2 


-18.5 > Mr 


0.025 <z < 0.05485 


74.6 < R< 162.6 


32550 (11341) 


3.78 


D3 


-19.0 > Mr 


0.025 <z < 0.06869 


74.6 < R< 203.0 


49571 (19270) 


4.18 


D4 


-19.5 > Mr 


0.025 <z < 0.08588 


74.6 < R < 252.9 


74688 (33039) 


4.58 


D5 


-20.0 > Mr 


0.025 <z < 0.10713 


74.6 < R < 314.0 


80479 (39333) 


5.56 



^Comoving distance in units of /i~^Mpc. 

''Number of early type galaxies 

•^Mean separation of galaxies in units of /i~^Mpc 



volvcd with a range of PSF sizes. For an observed galaxy 
image with a given size of the PSF we look for a best-fit 
Sersic model when the Sersic model is convolved with the 
PSF to find the true Sersic index and the true axis ratio. 
The fit is made at radii between 0.2 and l.Oiipet to avoid 
the central region, whose profile is strongly affected by 
seeing. When the best-fit Sersic models arc found in the 
g and i bands, the g — i color gradient and concentration 
index are calculated for the true Sersic models and their 
differences from those of convolved images are used to 
correct the color gradient of the observed image for the 
seeing effects. 

We classify galaxy morphologies into early (ellipticals 
and lenticulars) and late (spirals and irregulars) types 
using the automated galaxy morphological classification 
method given by Park & Choi (2005), who showed that 
incorporating the g — i color gradient into the classifica- 
tion parameter space allows us to successfully separate 
the two populations degenerated in it — r color and the 
concentration index helps us separate red disky spirals 
from early type galaxies (see Fig. 1 of Park & Choi 2005) . 
It should be noted that the "morphological" properties 
of galaxies can change when different bands are adopted. 
The difference in morphology in different bands is exactly 
the color distribution, which is reflected in our color gra- 
dient parameter. To define the classification boundaries 
in this multi-parameter space, they used a training set 
of 1982 SDSS galaxies whose morphological types are 
visually identified. The distributions of the early and 
late type galaxies in sample D2 are shown in Figures 6a 
and 66, below. The reader is referred to Park & Choi 
(2005) for further details of the methods and various 
tests. We have visually examined color images of a set of 
faint galaxies retrieved by the SDSS Image List, and con- 
firmed that the resulting morphological classification is 
highly successful with completeness and reliability reach- 
ing 90%. For the volume-limited samples D2 and CM 
we made an additional visual check to correct possible 
misclassifications of the automated scheme for blue early 
types (those below the straight line in Fig. 3a) and red 
late types (those redder than 2.4 in u — r color). The 
morphological types of 1.9% of galaxies, which are often 
blended or merged objects, were changed by this proce- 
dure. This morphology classification allows us to make 
very accurate measurements in this study compared with 
previous work. We discuss this further in section 3.2. 



3. RELATIONS AMONG PHYSICAL PARAMETERS 
OF GALAXIES 

We find that absolute magnitude and morphology are 

the most important parameters characterizing physical 
properties of galaxies, in the sense that other parameters 
show relatively little scatter once absolute magnitude and 
morphology are fixed. This is particularly true for early 
type galaxies. We begin this section by showing the de- 
pendence of various physical parameters of galaxies on 
their absolute magnitudes. In each case, we separately 
examine samples of early (E/SO) and late (S/Irr) mor- 
phological types. 

3.1. Variation with Absolute Magnitude 

Figure 3a shows the color-absolute magnitude diagram 

of galaxies in the CM and D2 samples. Galaxies in the 
CM sample are shown only for -17.5 > Mr > -18.45. 
The left panel shows the distribution of early types, and 
the right panel shows that of late types. The solid and 
dashed lines delineate the most probable u — r color of 
galaxies in the D2 and CM samples, respectively. (Here 
and throughout, the "most probable" parameter value 
within 0.4 absolute magnitude bin is the mode of the dis- 
tribution of that parameter.) The most probable color of 
the red sequence galaxies has a break at Mr w —19.6 and 
has slopes of about —0.07 at the bright side and about 
—0.15 at the faint side. We plot only those late type 
galaxies with axis ratios b/a > 0.6, to reduce the biases 
in absolute magnitude, color, color gradient, etc., due 
to internal extinction. Colors of late types have a much 
wider dispersion. Their most probable color rapidly be- 
comes bluer as galaxies become fainter, but is very in- 
sensitive to absolute magnitude fainter than Mr = —18, 
approaching u—r « 1.1 at the very faint end {Mr ~ —14, 
not shown). 

Baldry et al. (2004) reported a slope of about —0.08 for 
luminous red galaxies selected from the SDSS, which is 
consistent with our measurement. One significant dif- 
ference between Baldry et al.'s and our result is the 
color-magnitude relation of the late type (blue distribu- 
tion) galaxies at the bright side. In their results the 
slope of the relation becomes less steep long before the 
color-magnitude relation for late types touches that of 
early types as one, toward brightest magnitude. On the 
other hand, our result shows no such bending until it ap- 
proaches the relation for early types. The difference is 
probably caused by extinction of late type galaxies. 

Figure 3a shows that the color distribution of late type 
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Fig. 3. — Relations between the absolute magnitude and (a) u — r color, (fe) A{g — i) color gradient, and (c) Cin (inverse) concentration 
index of galaxies in the D2 (brighter than Mr = —18.5) and CM (fainter than —18.45) volume-limited subsamples of the SDSS. E/SO 
morphological types are shown in the left panels, and S/Irr types are shown in the right panels. Curves delineate the locations of the most 
probable parameter values as a function of absolute magnitude. The solid curve is for the D2 sample, and the dashed curve is for the CM 
sample. 



galaxies overlaps with early types at all magnitudes, and 
vice versa. The solid straight line in Figure 3a roughly 
divides the sequences of red and blue galaxies. How- 
ever, about 24% of spiral/irregular galaxies brighter than 
Mr = —18.5 are above this line when all inclinations are 
allowed. Among all galaxies above this line, 32% are 
actually late types. Therefore, an early type sample gen- 
erated by a simple cut in the absolute magnitude versus 
color space must take into account this huge contamina- 
tion. We also find a significant fraction of early types lo- 
cated outside the main red sequence at faint magnitudes. 
At absolute magnitudes between —18.5 and —19.0, about 
10% of early types are located below the straight line 
shown in Figure 3a. This outlying fraction is larger at 
fainter magnitudes. Some of these blue early type galax- 
ies may be actually spirals whose disks look too faint for 
both our automated and visual classification schemes to 
classify correctly, but the trend seems robust. 
Again, it should be noted that even though the bi- 



modality in galaxy color is mainly due to the bimodality 
in galaxy morphology, there is a significant overlap in 
color space between the two morphological classes. With 
few exceptions, very blue (u — r < 1.8) early type galax- 
ies show strong emission lines in their spectra, indicating 
active star formation (see Fig. 6c). The distributions of 
colors of early and late types at a given absolute mag- 
nitude are not Gaussian. The colors of early types are 
skewed to blue color, and those of late types are skewed 
to red color. The blue sequence appears less prominent 
when the data is contaminated by internal extinction ef- 
fects (see Fig. 12, below). 

Figure 36 shows distributions of galaxies in color gra- 
dient, A{g — i), versus absolute magnitude, Mr- We find 
that the majority of E/SO's {left panel) have a weak nega- 
tive color gradient; the envelope is slightly bluer than the 
core. Their most probable color gradient is surprisingly 
constant at A{g — i) = — 0.035 ± 0.007 at all magnitudes 
explored, with a weak trend for more positive gradient 
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toward faint and bright magnitudes. Interestingly, tlic 
E/SO's with Mr ~ Af, — 0.5 have the minimum color 
gradient. It should be noted that the color gradient of 
the brightest early types approaches zero, meaning that 
the brightest early types are reddest and most spatially 
homogeneous in color. Although the most probable color 
gradient has a nearly constant of absolute magnitude, the 
number of E/SO's with blue cores increases strongly at 
absolute magnitudes fainter than M^,. 

The color gradient of late type galaxies is more nega- 
tive (i.e., bluer envelope) than that of early types and is 
a slowly increasing function of absolute magnitude Faint 
late types are more gas rich and have more uniform inter- 
nal star formation activity, while bright late types have 
red central bulges and blue star forming disk structures. 

Figure 3c shows that the (inverse) concentration index 
Cin of the early types {left panel) is nearly independent 
of absolute magnitude when Mr < M*, but rises rather 
steeply at fainter magnitudes. The most probable con- 
centration index of the bright early types is close to that 
of a de Vaucouleurs profile (cjn = 0.29), even though 
it slightly increases at the brightest end. The brightest 
early types might be less concentrated because of more 
recent mergers, and the galaxies only slightly brighter 
than M*, with Mr ~ — 0.9 might be dynamically 
older (more relaxed). For late type galaxies the most 
probable value is close to the value an exponential disk 
(cin = 0.44) at the brightest end, albeit with a large dis- 
persion. 

The Petrosian radius in units oi kpc as a func- 
tion of absolute magnitude is shown in Figure 4o. The 
Petrosian radius is calculated by using the elliptical an- 
nul!, and is typically larger than the value in the SDSS 
Photometric database which uses the circular aimuli. 
Early type galaxies {left panel) exhibit a tight correla- 
tion between physical size and luminosity, and the rela- 
tion has a significant curvature. The slope is very steep 
{dlogRpet/dMr w -0.28 at Mr « -21.5) at the bright- 
est end. The late types have larger Rpet, and their size 
relation with absolute magnitude is less curved. The cur- 
vature of the relation for both early and late types means 
that fainter galaxies have lower surface brightness. As a 
guide the radius-magnitude relation in the case of con- 
stant surface brightness {L oc T,R^ where L is luminosity 
and E is surface brightness) is drawn in both panels. The 
slope of this relation to the bright end of the late types 
tells that the brightest spiral galaxies, with Mr < —21, 
have surface brightness approximately independent of lu- 
minosity. This relation is also true for the early types 
in the regime —21 < Mr < —20 fainter than for the 
late types. Desroches et al. (2006) also found that the 
radius-luminosity relation of elliptical galaxies is steeper 
at high luminosity. Shen et al. (2003) measured the size 
distributions of early and late type galaxies, divided by 
the concentration index or the Sersic index, and found 
that they are well fitted by a log normal distribution. 

Figure 46 shows the variation with absolute magni- 
tude of velocity dispersion, as measured by an automated 
spectroscopic pipeline called IDLSPEC2D version 5 (D. 
Schlegel et al. 2007, in preparation). Galaxy spectra 
of the SDSS are obtained by optical fibers with radius of 
iifib = 1-5"- The finite size of this aperture smooths the 
central velocity dispersion profile. To correct the central 
velocity dispersion for this smoothing effects we adopt 



a simple aperture correction, (Tcorr = crfib(8i?fib/-Ro)'' °^ 
(J0rgensen et al. 1995; Bernardi et al. 2003b), where Rq 

1 /2 

is the equivalent circular effective radius {b/a)^^yrdev 
where rdev is the seeing-corrected effective angular ra- 
dius along major axis of the galaxy from model fits to a 
de Vaucouleurs profile in the i band. Because the mea- 
sured velocity dispersion is systematically affected by the 
finite resolution of the spectrographs at small values, the 
calculation of the most probable velocity dispersion is re- 
stricted to galaxies with cr > 40 km s^"^. The errors in 
the measured velocity dispersions are about ~ 15% at 
70 km s^^ and ~ 30% at 50 km s~^. The dispersion in 
cr is typically Alogcr = 0.11, and the measurement er- 
ror is typically 11 km (see section 4.2). From these 
data we fit the slope, 7, defined by L « cr'*'. At the 
bright end, near absolute magnitude Mr = —22, we find 
7 = 4.6 ± 0.4, while the slope at magnitudes between 
Mr = -21 and -20 drops to 7 = 2.7 ± 0.2 (All five 
volume-limited samples are used for the measurements). 
For comparison, we plot two straight lines with slopes of 
7 = 4 {upper left, Fig. 46) and 3 {lower right, Fig. 46), 
the former being the slope of the Faber- Jackson relation. 
Recently, Desroches et al. (2006) also found that the 
slope 7 is steeper at high luminosity, which is consistent 
with our findings. The measured velocity dispersion of 
spiral galaxies brighter than Mr = -20.5 is 70% ± 2% 
of that of early types (or criate ~ Ceariy/\/2) at a given 
absolute magnitude. 

Because SDSS fiber spectra sample the light from only 
the central 1.5" radius region of galaxies, the measured 
spectral parameters are not fully representative of the 
physical properties of spiral galaxies, which tend to have 
large variation of star formation activity from the bulge 
to outer disk. Nevertheless, parameters based on fiber 
spectra, such as line widths and those from the principal 
component analysis, are often used to characterize the 
star formation activity, morphological types, etc. (e.g. 
Tanaka et al. 2004; Balogh et al. 2004b). Thus, for 
comparison, we show such results for our samples. 

The right panel of Figure 4c shows the cqiiivalent width 
of the Ha line of late types as a function of Mr- The 
equivalent width of the Ha line is often used as a measure 
of recent star formation activity (or nuclear activity). 
The majority of late types have large W{Ha). The plot 
shows that fainter galaxies tend to have larger W{Ha). 
This result should be considered carefully because the 
fiber spectra systematically miss the light from the outer 
disks of bright large galaxies, which might be active in 
star formation. The figure also shows that there is a 
class of late types that have weak Ha line emission. The 
placement in this diagram of some of these objects may 
be due to the finite fiber size, but some of them are gen- 
uine passive spirals (Couch et al. 1998; Poggianti et al. 
1999; Goto et al. 2003b; Yamauchi & Goto 2004). 

The left panel of Figure 4c shows W{Ha) of early type 
galaxies. We find that the number of early types with 
central star formation activity increases sharply at ab- 
solute magnitudes fainter than ~ M*. We find similar 
trends for m — r color and A{g — i) color gradient (see 
also Fig. 6, below). Early types with central star forma- 
tion activities have been studied and some of them are 
E -I- A galaxies (Goto et al. 2003c; Fukugita et al. 2004: 
Quintero et al. 2004; Goto 2005; Park & Choi 2005; 
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Fig. 4. — Relations between the absolute magnitude and (a) the Petrosian radius, (b) velocity dispersion, and (c) equivalent width of the 
Ha line of galaxies in the D2 (brighter than Mr = —18.5) and CM (fainter than —18.45) volume-limited subsamples of the SDSS survey. 
E/SO morphological types are shown in the left panels, and S/Irr types are shown in the right panels. Curves have the same meanings as 
in Fig. 3a — 3c. 



Yamauchi & Goto 2005). 

Figure 5 shows the isophotal axis ratio of early {left 
panel) and late (right panel) types as a function of ab- 
solute magnitude. We choose the isophotal axis-ratio 
rather than axis ratios from de Vaucouleurs or exponen- 
tial profile fits that are computed by the automated SDSS 
pipeline in our study because the isophotal position an- 
gles correspond most accurately to the true orientation 
of the major axis. The exponential profile fit is often bad 
particularly for barred spiral galaxies. Curves in Figure 
5 are the median axis ratios of galaxies in the D2 [solid 
line) and CM (dashed line) samples. We also draw the 
median curves in red and blue colors by using only galax- 
ies with b/a > 0.6. 

We clearly see that early type galaxies tend to be 
rounder at brighter magnitudes. Figure 5 demonstrates 
that the absolute magnitudes of inclined late types with 
b/a < 0.6 seriously suffer from dimming due to the inter- 
nal absorption, while the effects are ignorable for those 
with b/a > 0.6. See also Figure 12, below and discussion 



in section 4.3. 

3.2. Variation with Color 

Figure 6a shows the distribution of E/SO (left panel) 
and S/Irr (right panel) galaxies in the D2 sample in the 
u — r color versus lS.(g — i) color gradient space. Our 
morphological classification is based primarily on this pa- 
rameter space. Note that the sharp boundaries shown in 
Figure 6a are due to a choice of our morphological crite- 
ria in the u — r versus g — i color gradient space, which 
is determined to correlate well with visually classified 
morphologies. Scatter across the classification bound- 
aries is caused by the concentration index constraint on 
early types (see Park & Choi 2005 for details) and correc- 
tions made by visual inspection. Most early type galax- 
ies fall within a strong concentration near u — r — 2.8 
and /S.(g — i) — —0.04 in this plane. Because the color 
of the red sequence becomes bluer at fainter magnitudes 
(see Figure 3a) , the center of this concentration moves to 
bluer u — r when one examines fainter early types. There 
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Fig. 5. — Isophotal axis ratios of the early (left) and late (right) type galaxies in sample D2 (solid magenta line) and CM (dashed green 
line). Curves are the median values. For late types we added the median curves, in red and blue colors, derived from only those with 
b/a > 0.6. 



is also a trail of early types towards bluer colors and 
more positive color gradients. The classification bound- 
ary separating this blue early types from late types has 
been determined from a training set that showed such 
a trend (Park & Choi 2005). The majority of galaxies 
in this early type trail are fainter than (6Zue points), 
often show emission lines (see early types in Fig. 6c), 
and live mainly at intermediate and low density environ- 
ments (see early types in Fig. 11a of Paper II). These 
blue early types amount to 10% of early types, which af- 
fects the morphology fraction as a function of luminosity, 
the type-specific luminosity functions, and so on. In this 
manner, the color gradient and concentration constraints 
in our morphological classification reduce the misclassifi- 
cation of red late type as early type, and blue early type 
as late type. This is one of the major improvements we 
made in this study compared with previous works. Late 
type galaxies brighter than Af, are redder and have more 
negative color gradients than fainter late types. Figure 
6b shows galaxies in the u — r versus the (inverse) con- 
centration index space. The brightest early types have 
li — r sa 2.9 and Qn ~ 0.3. There is a weak tendency for 
the redder (2.6 <u — r< 3.0) bright early types to have 
smaller Qn than the bluer (2.3 < u—r < 2.6) fainter ones. 
Late types are loosely concentrated at u — r « 1.5 and 
Cin = 0.49. They have a broad overlap with the region 
occupied by early type galaxies. 

Figure 6c shows the relations between u — r color and 
the equivalent width of Ha line, W{Ha). Most red 
early type galaxies have weak emission, W{Ha) < 3. 
Because galaxies in the blue trail of early types usually 
have blue star-forming centers, W{Ha) in the light from 
their centers has a tight correlation with their color. The 
right panel of Figure 6c shows that the Ha emission line 
strength of the centers of late type galaxies (plotted as 
log[W^(i?a)]) is close to being linearly proportional to the 
overall u — r color. This is particularly true for blue late 
types. This proportionality is as expected if Ha emission 
is proportional to it-band luminosity. It is interesting to 
observe that the sequence formed by the blue early types 
is far from a straight line and has higher W{Ha) than 
the sequence of the blue late types. 

3.3. Velocity Dispersion of Early Types 



Stellar velocity dispersion is a dynamical indicator of 
the mass of galaxies. Figure 7a shows that there is a 
good correlation between velocity dispersion and concen- 
tration index, and thus between mass and surface bright- 
ness profile. Here it should be noted that the D2 sample 
starts to be incomplete at < 150 km s^^ due to the 
absolute magnitude cut and the significant dispersion of 
cr at a given absolute magnitude (see Fig. 46). 

The relation becomes tighter for high-cr early type 
galaxies and approaches the de Vaucouleurs profile value 
of Cin = 0.29 at the highest end. The size of galaxies is 
well correlated with velocity dispersion when a > 100 
km s~^, but the dependence become weaker at lower 
(T (Figure 7b). Note that the number of galaxies with 
Rpet ^ 3/i~^kpc is underestimated in the D2 sample. 
The error in a is typically 11 km (see below), which 
somewhat broadens the size versus a relation. These 
systematic effects should be taken into account in inter- 
preting Figure 7b. 

Figure 7c indicates that there are progressively more 
centrally star-forming early type galaxies as the veloc- 
ity dispersion decreases, consistent with the relation be- 
tween W{Ha) and AU (Fig. 4c). 

4. COLLECTIVE PHYSICAL PARAMETERS 

4.1. Luminosity Function 

The upper panel of Figure 8 shows the luminosity func- 
tion (LF) of galaxies in the D3 sample. Filled circles are 
the LF measured by using both early and late type galax- 
ies. Open circles are the LF of the E/SO morphological 
types, and open squares are that of the S/Irr types. Sym- 
bols are calculated by binning galaxies in each magnitude 
bin. Uncertainty limits are estimated from 16 subsets of 
each sample. Curves are the best-fit Schechter functions 
of the following form 

(j){L)dL = (j)., {L/L^)°'exp{-L/L^)dL/L^. (3) 

We use the MINUIT package of the CERN program 
library^ to determine the parameters in the Schechter 
function by the maximum likelihood method, which is 
applied to individual galaxies as described in Sandage 

http:/ /http:/ /wwwasdoc. web. cern.ch/wwwasdoc/minuit/minmain. html 
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Fig. 6. — Relations between the u — r color and (a) g — i color gradient, (b) (inverse) concentration index, and (c) the equivalent width 
of the Ha line. The most likely Ha equivalent width is drawn for the late types (lower right). Blue points are fainter galaxies with 
Mr > —20.0, and red points are those with Mr < —20.0. 



et al. (1979). The Schechter function fits the data ex- 
tremely well except for the highest luminosity bins. To 
reduce the bias in LF due to the internal absorption, 
we have included only those late types with axis ratios 
b/a > 0.6 in our analyses. To estimate the number den- 
sity (or 0*) of all spiral galaxies, we need to know the 
fraction of late type galaxies with apparent axis ratios 
greater than 0.6. We adopt the fraction 0.505 calculated 
from the CM sample (see Fig. 135, below), whose spiral 
membership is considered to be the least affected by the 
internal absorption. We use this fraction to scale up the 
late type galaxy LF estimated by using only those with 
b/a > 0.6. The results of the Schechter function fits are 
summarized in Table 2. 

Since our volume-limited samples from the SDSS do 
not extend to very faint magnitudes, there is a signifi- 
cant correlation between measured M, and a. We choose 
the D3 sample to estimate a representative LF of the 
SDSS galaxies as a compromise; the uncertainty in 
is minimized for D4, while the uncertainty in a is min- 
imized for D2. However, the D3 sample covers the ab- 



solute magnitude space brighter than —19.0, only about 
1.2 mag fainter than M*, and the correlation between 
the estimated parameters M^, and a makes their physi- 
cal meanings less clear. The LF of the deeper D4 sample, 
which includes the Sloan Great Wall (Gott et al. 2005), 
shows the highest density of early type galaxies among 
all five volume-limited samples, and thus its LF seems 
less representative. The D5 sample does not extend faint 
enough to determine the a parameter. The existence of 
the large scale structure affects both normalization and 
shape of LF (e.g. Croton et al. 2005). Therefore, the 
trends in and a are affected by the fact that nearby 
universe is a relatively underdense region and the shal- 
lower samples are more dominated by the relatively faint 
galaxies which prefers underdense regions (see Fig. 10 
of Paper II for further details). When all spiral galaxies 
are used (i.e. no internal absorption effect correction), 
the LF of all galaxies in the D3 sample is well fitted by 
M, = -20.24 ± 0.03 and a = -0.90 ± 0.03. For compar- 
ison, Blanton et al. (2003a) report A/, = -20.44 ± 0.01, 
and a = —1.05 ±0.01 for an apparent magnitude-limited 



Properties of SDSS Galaxies 



11 



TABLE 2 

Total and Type-specific Schechter Function Parameters 



Name 


Absolute Magnitude 


M, - 51ogio/i 


a 


4>* 

10-2/l3Mpc-3 


All Types: 










Dl 


-18.0 > Mr 


-20.31 ± 0.04 


-0.918 it 0.027 


1.64 


D2 


-18.5 > Mr 


-20.22 ± 0.03 


-0.843 it 0.026 


1.72 


D3 


— 19.0 > Mr 


—20.22 ± 0.03 


—0.807 it 0.028 


1.85 


iJ4 


— ly.o > IVIr 


on oc _L n no 
— zu.zo it U.Uz 


n Qii'y J- A nun 


O A1 


Do 


— 2U.U > Mr 


OA QO _l_ A no 


n OAK _1_ A AQA 


1.87 


Early Types: 










Dl 


-18.0 > Mr 


-20.17 it 0.05 


-0.494 it 0.047 


0.73 


D2 


— 18.5 > Mr 


—20.16 it 0.04 


—0.473 it 0.044 


0.71 


D3 


-19.0 > Mr 


-20.23 ± 0.04 


-0.527 it 0.043 


0.71 


D4 


-19.5 > Mr 


-20.38 it 0.03 


-0.719 it 0.042 


0.75 


D5 


-20.0 > Mr 


-20.49 it 0.03 


-0.870 it 0.052 


0.66 


Late Types: 










Dl 


-18.0 > Mr 


-20.30 it 0.05 


-1.070 it 0.033 


0.98 


D2 


-18.5 > Mr 


-20.18 it 0.04 


-0.965 it 0.034 


1.10 


D3 


-19.0 > Mr 


-20.12 it 0.03 


-0.897 it 0.037 


1.23 


D4 


-19.5 > Mr 


-20.09 it 0.03 


-0.832 it 0.043 


1.44 


D5 


-20.0 > Mr 


-20.04 it 0.03 


-0.666 it 0.061 


1.37 



sample called SDSS LSS Sample 10 (similar to SDSS 
DRl). If we adopt the faint end slope a from the shal- 
low sample Dl and from the deepest sample D5, our 
estimates arc a = -0.99 ±0.03 and Af* = -20.38 ±0.03, 
which are vc;ry close to Blanton et al.'s. 

When the internal absorption effect in the sample is 
reduced, the LF of spirals changes, particularly in the 
a parameter. The LF of early types in sample D3 has 
= -20.23 ± 0.04 and a = -0.53 ± 0.04 while that 
of late types has a fainter characteristic magnitude of 
M* = —20.12 ± 0.03 and a steeper faint-end slope of 
a = — 0.90±0.04. If all spiral galaxies are used, the LF of 
late types has M* = -20.09 ±0.04 and a = -1.01±0.04. 
This demonstrates that the faint-end slope a of late types 
would have been misleadingly measured steeper if the 
inclined late type galaxies had not been excluded. The 
blue galaxies classified as early types by our morpho- 
logical classifier make an increasing contribution to the 
LF at faint magnitudes. We suspect that some of these 
objects are actually bulge-dominated spirals whose disks 
are lost to the background sky. However, the absolute 
contribution of blue early type galaxies to the LF is still 
small at the magnitudes studied, as seen in Figure 3a, 
and the type-specific LFs at magnitudes brighter than 
Mr — 18 seems well determined by our samples. Baldry 
et al. (2004) used u — r color-selected early (red) and 
late (blue) type galaxy samples to measure their lumi- 
nosity functions. The estimated M* parameter differed 
by 0.21 mag and the a-parameter was measured to be 
-0.83 (early) and -1.18 (late). The difference in Af* be- 
tween early and late types is quite consistent with our 
result, but they obtained steeper a's. 

The bottom panel of Figure 8 shows the fraction of 
E/SO galaxies in all five volume-limited samples as a 
function of r-band absolute magnitude. Again, only late 
types with b/a > 0.6 are used, and their number frac- 
tion of 0.505 is used to infer the total number of spiral 
galaxies. There is a surprisingly consistent monotonic 
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Flc;. 7. — Relation between the velocity dispersion of early type 
galaxies in sample D2 and (a) (inverse) concentration index, (6) 
the Petrosian radius, and (c) the equivalent width of the Ha line. 
Larger dots are galaxies with Mr < —20.5, and smaller dots are 
those with Mr > -20.5. 
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Fig. 8. — Top: Luminosity function (LF) of galaxies in the D3 
sample. Filled circles are for all galaxies, open circles for the E/SO 
morphological types, and open squares for the S /Irr types. Symbols 
are the LF calculated by binning galaxies in each magnitude bin. 
Curves are the best-fit Schechter functions. Bottom: The elliptical 
galaxy fraction in magnitude bins measured from the five volume- 
limited samples (Dl - D5). 

relation between luminosity and early type galaxy frae- 
tion. The critical magnitude at which the early type frac- 
tion exceeds 50% is = -21.29 ± 0.04. Even though 
the shape of LF itself shows significant fluctuations due 
to large scale structures when measured in different re- 
gions of the universe, the morphological fraction as a 
function of luminosity is relatively less sensitive to large 
scale structures and thus seems to be more universal. 
However, it is well-known that the early type fraction 
is a monotonically increasing function of local density 
(Dressier 1980; Tanaka ct al. 2005; Goto et al. 2003a; 
Paper II). Weinmann ct al. (2006) also measured the 
morphology fraction as a function of Mr and found basi- 
cally the same dependence, even though the sample they 
used is smaller. 

4.2. Velocity Dispersion Function 

The velocity dispersion function (VDF) of early type 
galaxies can be used to derive gravitational lens statis- 
tics and to estimate cosmological parameters. To esti- 
mate the velocity dispersion function we need to have a 
sample of galaxy velocity dispersion data that is com- 
plete over a wide range. Figure Ah indicates that sample 
D2, for example, starts to be incomplete at ct < 150 
km s~^ due to the absolute magnitude cut and the sig- 
nificant dispersion of ct at a given absolute magnitude. 
To lower the completeness limit down to ct w 70 km 
we have used a series of our volume-limited sam- 
ples an additional very faint volume- limited sample with 
— 16.8 > Mr > —17.5, which is added because even the 
CM sample is not faint enough to give the CT-distribution 
complete down to ct = 70 km s^^. The differences in 
volumes of these volume-limited samples are taken into 
account when we estimate the velocity dispersion func- 
tion from a combination of seven samples. 
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Fig. 9. — Velocity dispersion functions (VDFs; circles) measured 
from our samples with galaxies brighter than Mr = —16.8. The 
solid curve is our estimate of the true VDF best fit by a modified 
Schechter function. The long-dashed curve is the best fit to the 
raw measurement. The dotted curve is the fit given by Sheth et al. 
(2003). 



o 



10* 



5000 



2 



T" 



T" 



nJ 



Jl 



T" 



IDLSPEC2D 
IDLSPEC2D repeat 
measurements 

SPECTROID repeat 
measurements 



T" 



20 40 
Error (ct) (km s~') 



60 



Fig. 10. — Histogram showing the distributions of the standard 
deviations of a of the galaxies measured by two different SDSS 
pipelines IDLSPEC2D {thin and thick solid lines) and SPECTROID 
{dotted line). The thick solid and dotted lines are obtained from 
the repeat measurements (plate 406 and 412) and the thin solid 
line is from the formal error given by the IDLSPEC2D. 

The circles in Figure 9 show the VDF of early type 
galaxies measured in this way. The curve is the fit by 
a modified Schechter function of the form 

/3 da 



(4) 



'l](a//3) CT 

where ct* is the characteristic velocity dispersion, a is the 
low-velocity power-low index, and /? is the high-velocity 
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exponential cutoff index (Mitchell et al. 2005). The VDF 
parameters obtained by the MINUIT package are 

((?!)*, C7*, a, /3) = (8.0 X IQ-^/i^Mpc"^, 161 ± Skms'S 
2.32 ±0.10, 2.67 ±0.07). 

When we made the fit with the above parameters, we 
assumed that the measurement error in cr was approxi- 
mated by 0.0840- km s^^ (sec below). The dashed curve 
in Figure 9 is the modified Schechter function best fit to 
the raw data points {open circles) which are broadened at 
high a's due to the measurement error. The solid curve is 
the deconvolved VDF whose parameters are given above. 

It should be noted that there exists very strong corre- 
lations among the fitting parameters, and the results of 
the fitting depend sensitively on the range of the velocity 
dispersion used. Our estimate of the VDF of early type 
galaxies is quite different from that of Sheth et al. (2003) 
and Mitchell et al. (2005). The dotted curve in Figure 9 
is the VDF best fit given by Sheth et al. (2003). Mitchell 
et al. (2005) presented the same function with the nor- 
malization lowered by 30%. The difference at small a 
seems due to the fact that some of morphologically early 
type galaxies are discarded based on spectrum in Sheth 
et al. (cf. Fig. 12 of Mitchell et al. 2005). It is not clear 
what makes the difference at large a, but it could be due 
to the measurement error larger than that adopted here. 

In order to check whether or not our VDF is seriously 
affected by the measurement error in a we have made 
the following test. We first adopt Sheth et al. (2003)'s 
VDF as the true one. We then look for galaxies hav- 
ing 4 or 5 repeat measurements of a in southern stripe 
82 (plates 406 and 412) to estimate the error. The thick 
histogram in Figure 10 shows the distribution of the stan- 
dard deviations of cr of the galaxies with repeat measure- 
ments. The median size of error is 11 km . One can 
note that the error estimated from the repeat measure- 
ments by the IDLSPEC2D is larger than the formal error 
[thin histogram) given by the IDLSPEC2D, but is compa- 
rable with that (dotted histogram) obtained from the re- 
peat measurements by the official spectroscopic pipeline, 
SPECTROID (written by M. SubbaRao, M. Bernardi and 
J. Frieman). If the measurement error were larger than 
20 km s""'^, the dispersion of the relation between Mr 
and a shown in Figure 4& would be much larger. We 
then calculate the VDF affected by the error in cr. The 
dotted curve in Figure 11 is Sheth et al.'s function, and 
the short-dashed curve is the VDF deformed due to the 
Gaussian-distributed measurement error of 0.084cr km 
s"^, which reflects the fact that the error in a estimated 
from the repeat measurements increases as a increases. 
It demonstrates that our VDF cannot be realized from 
Sheth et al.'s function by the convolution effects of the 
measurement error 

4.3. Axis Ratio Distribution 

The scatter plot of galaxies in Figure 5 indicates that 
the mean axis ratio of bright late type galaxies tend to 
be high. This is due to the internal absorption making 
inclined galaxies appear fainter. The bottom panel of 
Figure 12 shows the late type galaxies of sample D3 in 
the Mr versus u — r color space as in the right panel of 
Figure 3a. Galaxies with b/a > 0.8 are shown as blue 
dots, and those with b/a < 0.4 are shown as red dots. 
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Fig. 11. — Velocity dispersion functions (VDFs) demonstrating 
the effect of the measurement errors in cr on the VDF. The dotted 
curve is Sheth et al.(2003)'s VDF, and the short-dashed curve is 
the VDF deformed due to the Gaussian-distributed measurement 
error of 0.084cr km s^^. The long-dashed curve is the best fit to 
the VDF measured in this work. 

Their most likely color curves are also drawn. We find 
that the group of highly inclined galaxies is significantly 
shifted with respect to the almost face-on ones toward 
fainter magnitudes and redder colors. We also find that 
the shapes of the two distributions are different, which 
can be explained if the galaxies with intermediate mag- 
nitude {Mr ^ —20.5) and color {u — r ^ 1.6) sufi^er from 
more dimming and reddening due to internal extinction. 
The color of very red or very blue galaxies seems to be 
less affected by the internal extinction. A detailed study 
of internal extinction is needed for the SDSS galaxies. 
The early types are not much affected. If one does not 
take into account this inclination effect on magnitude, 
the measured LF of late type galaxies and the morphol- 
ogy fraction will be seriously in error. In analyses where 
these systematics are important, we have used late types 
with b/a > 0.6. Figure 5 demonstrates that the inclina- 
tion efiFects on magnitude is small for this set of late type 
galaxies. 

Because of the inclination effects on magnitude, our 
volume-limited samples are not complete for late type 
galaxies near the absolute magnitude cuts. The lower 
panel of Figure 13 shows distributions of isophotal axis 
ratio of late type galaxies in samples CM (solid line with 
filled circles), D2 {short-dashed line with open circles), 
and D4 {long-dashed line with open squares). There ap- 
pears to be no late type galaxy with nearly zero axis 
ratio because of the intrinsic thickness of disks of late 
type galaxies. The axis ratio of 0.2 due to the intrin- 
sic thickness is consistent with previous works (Sheth 
et al. 2003). Almost face-on galaxies are also lacking, 
because disks of late type galaxies do not appear ex- 
actly round even when they are face-on (Kuehn & Ry- 
den 2005). The distributions are significantly different 
from one another because of the internal absorption ef- 
fects, lacking highly inclined galaxies in the samples with 
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Fig. 12. — Distributions of (a) early and {b) late type galaxies in the color magnitude diagram. Galaxies with b/a > 0.8 arc shown 
as blue dots, and those with b/a < 0.4 are shown as red dots. For late types the most likely color-magnitude relations are drawn for the 
different axis ratio sets. 



bright cuts. We assume that the magnitude hmit of CM 
sample (M^^cut = —17.5) is faint enough that its b/a 
distribution is close to the intrinsic distribution. The 
fraction of galaxies with b/a > 0.6 in the CM sample 
is 0.505. The upper panel of Figure 13 shows the dis- 
tribution of projected axis ratio of early type galaxies. 
Distributions obtained from three samples with differ- 
ent magnitude cuts match one another very well. This 
reflects the fact that internal absorption in early type 
galaxies is negligible. 

Alam & Ryden (2002) measured the axis ratio distri- 
butions for red and blue populations in the SDSS EDR 
sample. The distribution functions for red or blue galax- 
ies well-fitted by a de Vaucouleurs profile match well 
with our result for early types. Sheth et al. (2003) also 
measured the distribution of axis ratio of early and late 
type galaxies in the SDSS. Their axis ratio distribution 
of early types agrees well with ours except at the b/a = I 



bin. 

5. SUMMARY 

We have elaborately investigated the relations among 
various internal (e.g. color, luminosity, morphology, star 
formation rate, velocity dispersion, size, radial color gra- 
dient, and axis ratio) and collective properties (e.g. lu- 
minosity function, velocity dispersion distribution func- 
tion, and axis ratio distribution function). We will study 
other collective properties such as the correlation func- 
tion, power spectrum, topology, and peculiar velocity 
field, in separate works. 

There are a number of important improvements we 
made in our measurement: 

1. To extend the absolute magnitude range, we have 
added the redshifts of the bright galaxies with 
r < 14.5 to the SDSS-NYU-VAGC. Six volume- 
limited samples with faint limits from Mr = —17.5 
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Fig. 13. — Distributions of the apparent isophotal axis ratio of 
(a) early and (6) late type galaxies. 

to —20.0, are used for more efficient use of the flux- 
limited SDSS sample. 

2. We have divided the samples into subsamples of 
early (E/SO) and late (S/lrr) morphology types by 
using our accurate morphology classifier, working 
in the three-dimensional parameter space of m — r 
color, g — i color gradient, and concentration index. 
The morphology classifier is able to separate the 
blue early types from the blue late types, and also 
to distinguish the spirals from the red E/SO types 
because it does not depend only on color. 

3. The late type galaxies with isophotal axis ratio of 
less than 0.6 have been excluded because they suf- 
fer severely from dimming and reddening due to 
internal extinction and then may cause biases in 
luminosity, luminosity function, color, color gradi- 
ent, and so on. This correction makes the blue se- 
quence of the late type galaxies bluer and brighter 
compared to that obtained from all late type galax- 
ies. The blue sequence also become better-defined 
with smaller width. 

We have found that absolute magnitude and morphol- 
ogy are the most important parameters in characterizing 
physical properties of galaxies. In other words, other pa- 
rameters show relatively small scatter once absolute mag- 
nitude and morphology are fixed for early type galaxies, 
in particular. One very interesting fact we have noticed 
in our work is that many physical parameters of galaxies 
manifest clifFercnt behaviors across the absolute magni- 
tude of about M* ± 1. For example, the red sequence 
of early types changes the slope at Mr « —19.6 in the 
color-magnitude space. Also, at a magnitude fainter than 
Mr ~ —20, the number of blue star-forming early types 
increases significantly, and the surface brightness profile 



as parameterized by Cin becomes less centrally concen- 
trated. The passive spirals with vanishing Ha emission 
seem to have magnitudes brighter than ~ —19.5. 

At fixed morphology and luminosity, we find that 
bright {Mr < —20) early type galaxies show very small 
dispersions in color, color gradient, concentration, size, 
and velocity dispersion. These dispersions increase at 
fainter magnitudes, where the fraction of blue star- 
forming early types increases. Concentration indices of 
early types are well-correlated with velocity dispersion, 
but are insensitive to luminosity and color for bright 
galaxies, in particular. The slope of the Faber- Jackson 
relation (L cx a'^) continuously changes from 7 = 4.6±0.4 
to 2.7 ± 0.2 when luminosity changes from Mr = —22 to 
—20. The size of early types is well-correlated with stellar 
velocity dispersion, cr, when a > 100 km s~^. Late type 
galaxies show wider dispersions in all physical parame- 
ters compared to early types at the same luminosity. We 
find that passive spiral galaxies are well-separated from 
star-forming late type galaxies at Ha equivalent width 
of about 4. 

We note that the estimated LF shows significant fluc- 
tuations due to large scale structures when measured in 
different regions of the universe. On the other hand, the 
morphological fraction as a function of luminosity is rel- 
atively less sensitive to large scale structures and thus 
seems to be more universal. In Paper II, it is shown 
that the morphology fraction is a monotonic function of 
local density at a fixed luminosity. Since luminosity in 
turn depends on local density, this result indicates that 
the probability for a galaxy to be born as a particular 
morphological type is basically determined when its lu- 
minosity or mass is given. To this extent, morphology is 
determined by galaxy mass. The question here is what 
determines the morphology of a galaxy having a partic- 
ular mass while keeping the type fraction corresponding 
to its mass scale. Further investigations are required to 
answer this question. 
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